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У роботі розглянуто математичні моделі біосенсорних та імуносенсорних динамічних систем у вигляді 
диференціальних рівнянь у частинних похідних. Запропоновано моделі для оптимізації розробки біосенсорів, 
модель біосенсора в циліндричних координатах, на основі використання кінетики Міхаеліса-Ментена та 
рівнянь реакції-дифузії. Розроблено модель імуносенсора у вигляді решітчастих деференціальних рівнянь із 
запізненням. Побудова моделі грунтується на ряді біологічних припущень щодо взаємодії колоній антигенів та 
антитіл, а також дифузії антигенів. Для опису дискретних у просторі колоній, локалізованих у відповідних 
пікселях, використовується апарат решітчастих диференціальних рівнянь. 
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Вступ 

Протягом останніх років спостеріга- 
ємо стрімкий розвиток науки і техніки, що 
зумовило появу нових методів детекції. 
Тому в науці та промисловості зростає 
інтерес до біосенсорів. Біосенсори є 
альтернативою до загальновикористову- 
ваних методів вимірювання, які характе- 
ризуються поганою вибірковістю, високою 
вартістю, поганою стійкістю, низьким від- 
гуком і переважно можуть використову- 
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ватися лише високодосвідченим персона- 
лом. Біосенсори - це аналітичні прилади, 
які перетворюють процеси біохімічних ре- 
акцій на фізичну величину, яку можна ви- 
міряти. У своїй конструкції вони викори- 
стовують біологічні матеріали, які надають 
високу вибірковість, селективність, точ- 
ність, дають змогу здійснювати швидкі і 
прості вимірювання |. 

Біосенсори характеризуються висо- 
кою ефективністю і широко використову- 
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ються у харчовій промисловості |2|, при 
захисті навколишнього середовища |З, в 
оборонній промисловості (4|, але найчас- 
тіше - у медицині |35-8|, як інструмент для 
постановки діагнозів. У цілому, сімейство 
біосенсорів ділиться на дві частини. Перша 
пов'язана з рівнем наближення рецептора 
до біологічного матеріалу, який викорис- 
товується в його будові. Рецепторами мо- 
жуть бути ензим, протеїн, порферін, анти- 
ген або антитіло. Друга частина біосен- 
сорів обмежена до шару провідника, де 
біологічний ефект перетворюється на 
вимірювальний сигнал, який може бути 
електрохімічний, імпедансний, амперомет- 
ричний, оптичний та ін. 

Постановка проблеми 

Математичні моделі біосенсорів та 
імуносенсорів повинні враховувати про- 
сторово-часові властивості пристроїв, у 
яких використовується детектор. Відносно 
просторової організації досліджувана мо- 
дель повинна грунтуватися на певній 
дискретній структурі, яка буде врахо- 
вувати взаємодію пікселів імуносенсора. З 
точки зору часових змін, моделі повинні 
описувати процеси, відомі як популяційна 
динаміка. Саме тому проблему становить 
розробка математичних моделей біосен- 
сорів та імуносенсорів, які б враховували 
біологічні припущення щодо основних 
компонент пристрою - симетричну геомет- 
рію біосенсора та однорідний розподіл 
іммобілізованого ферменту та реакції 
зв'язування у ферментному шарі, кінетичні 
властивості бактеріальних  ферментатив- 
них реакцій, кінетику Міхаеліса-Ментена, 
колоній антигенів і антитіл, локалізованих 
у пікселях, дифузії колонії антигенів між 
пікселями та ін. 

Аналіз 
публікацій 

Останнім часом велика увага дослід- 
ників привернута питанням розробки та 
використання біосенсорів та імуносенсорів 
П-10). У роботі (|5| наведено грунтовний 
огляд теоретичних основ проектування 
біосенсорів. Прикладні аспекти використання 
імуносенсорів обговорюються в |З). 

Важливе значення у функціонуванні 


останніх досліджень |і 
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біосенсорів має фізичне явище флуорис- 
ценції, про що йде мова в |3, 4, 10). У (13| 
сформульовано основні задачі, пов'язані з 
дослідженням стійкості у біосенсорах. У 
роботах |, 11, 12| йдеться про проек- 
тування сенсорів, в основу роботи яких 
покладено перебіг хімічних реакцій. У 
роботі |12| для такого роду сенсорів було 
запропоновано математичне моделювання 
в класі решітчастих диференціальних 
рівнянь. У даній роботі такий підхід буде 
використано для імуносенсорів. 

У П4) викладено основні результати 
щодо використань рівнянь популяційної 
динаміки. У роботі |11| розроблена і дос- 
ліджується спрощена модель імунної сис- 
теми за допомогою диференціальних рів- 
нянь із запізненням, що буде використана 
в даній роботі. 

Мета дослідження 

Запропонувати моделі біосенсорів та 
імуносенсорів у вигляді диференціальних 
рівнянь у частинних похідних. 

1. Структура  біосенсорів та 
імуносенсорів 

Біосенсор - це аналітичний прилад, 
що містить у своєму складі біологічний 
чутливий елемент (фермент, антитіло, 
ДНК, клітинні органели, клітини чи шма- 
точки тканин), поєднаний з перетворю- 
вачем (електрохімічним, оптичним, кало- 
риметричним, акустичним). Вимірювання 
концентрації мішені аналізатор виконує 
кількісним перетворенням параметрів ре- 
акції у кількісний електричний чи оптич- 
ний сигнал. 

Серед великого сімейства біосенсорів 
імуносенсори є типовими сенсорами, що 
містять шар рецептора, який чутливий 1 
селективний, включаючи  імобілізований 
біологічний елемент, наприклад, антитіло, 
антиген або хаптен, які є імунологічними 
рецепторами для молекул, які дослід- 
жуються. В імунсорі (імуносенсорі) відбу- 
вається реакція, яка грунтується на взаємо- 
дії між антитілом і антигеном або малень- 
кими молекулами (хаптенами). Антитіла 
часто називаються імуноглобулінами тому, 
що вони є протеїнами, які пов'язані з 
імунною системою. 
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Імуноглобуліни | використовуються 
імунною системою для ідентифікації та 
нейтралізації чужорідних об'єктів. Вони 
використовують властивості зв'язування 
антигенів. Антигени і антитіла можуть 
використовуватися в шарі рецептора у біо- 
сенсорах. Зменшення властивостей, які 
пов'язані з антитілами під час процесу 
імубілізації антигена, використовуються в 
конструкції шару рецептора, де антитіла 
відіграють функцію  аналітів (молекул 
предметного детектування) |З9). 

Молекули, які пов'язані з детекцією, 
забезпечують зв'язування антитіл з анти- 
генами, утворюючи складні конструкції. 
При цьому між антигенами і антитілами 
утворюються дуже сильні зв'язки з кон- 
стантою зв'язування Ка -10 -10"" (10). 

2. Моделі біосенсорів у вигляді 
диференціальних рівнянь в частинних 
похідних 

21. Модель для 
розробки біосенсорів 

Оптимізаційні методи проектування 
(15) біосенсорів вимагають математичного 
моделювання. Припускаючи симетричну 
геометрію  біосенсора | і | однорідний 
розподіл іммобілізованого ферменту та 
реакції зв'язування у ферментному шарі, з 
одновимірною в просторі дифузією, яка 
описана другим законом Фіка, наведемо 
рівняння реакції такого типу дифузії 
(О«а«а, 1»0). 


оптимізації 


дЕ д'Е 
З СДН р ЗПС ГО ГУ ОА п 

ді тн де СН г А 
Бу 2 а! абз 
212 ОРЕ 

оц) о о КЕ СЯ 

ді яю да Ре р (2) 
ч ЖЕ а Я 292 
7 ОР о 
кою во іх КО 
КР, Як Р28» 
05, М 978, 
сли он о 


нК,5, Р, - Кк, Р5, 
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др, д'Р. 

а ка по З (5) 

кВ, Кк В5, 

ОР, д'Р. 

-асзар 252 Е - 

ОЇ Р, да? ня КУ сарбз ; (6) 

-кР кр, 

дс Же 

тю Р; зона ОЮ (7) 
де с 1 1 означають простір і час, 
відповідно, ЩЕ, (2, г) і ЕЕ, (2, г) і 


концентрації окисленої (Е, ) і зменшеної 


(Ем) глюкозодегідрогенази, відповідно, 
5 (2,г) і 5. (51) - концентрації субстратів 
(фериціаніду та окисленого медіатора), 
Р(21) і Рог) - концентрації продуктів 
реакції  (фериціаніду | і | зменшеного 
медіатора), С(2,1) концентрація глю- 
кози, ДВ. ,, Де, , В,» Др, Р. - відповідні 
коефіцієнти дифузії в ферментному шарі, а 
а, - товщина ферментного шару, і- 1,2. 
За межами шару ферменту відсутня 


ферментативна реакція 124, Кінетика 


неферментативних реакцій так само як 
транспорт масопередачі через мембрану 


(4 «««ача,) ; у зовнішньому шарі 
(4, ча, «с«а ча, ча,) описується 


аналогічними системами рівнянь реакції- 
дифузії 16). 

Спочатку фермент 5, рівномірно 
поширюється | у  ферментному шарі 
О«сх«а,. Біосенсор починає працювати 
при (1-0), коли з'являється субстрат 5, 
та глюкоза на зовнішній межі дифузійного 
шару (с-4, -4., 4). 

Рівняння керування (1-7) є типовими 
для класу біосенсорів, які використо- 
вуються для визначення субстратів 1171. 

2.2. Модель біосенсора в 
циліндричних координатах 

У роботі |13| розглянуто матема- 
тичну модель біосенсора, що описує кіне- 
тичні властивості бактеріальних фермента- 
тивних реакцій у відповідь на токсичні 
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хімічні речовини та отриману електро- 
активну дифузію молекули в мініатюрній 
електрохімічній комірці електроду. Мо- 
дель характеризує генерований електрич- 
ний струм як функцію бактерій та 
концентрацій токсинів, електрохімічних 
розмірів комірок та електроду. 

Для того щоб розрахувати поточний 
струм як функцію часу, зімітовано дифу- 
зію молекул електрохімічного активного 
продукту ферментативної реакції з кон- 
центрацією С, (г, 2,1). Ці молекули окис- 
люються на робочій поверхні електрода та 
генерують струм. У моделі зроблено при- 
пущення, що струм повністю контро- 
люється масою переносу та кінетикою 
реакцій у розчині на електроді. 

Оскільки електрохімічна комірка має 
циліндричну симетрію, то можна написати 
рівняння дифузії в циліндричних коорди- 
натах, відповідно до другого закону Фіка: 


ОС (г, 2,1) -р еекуво рт 
аинрчанноії 2/4 Порча 


ОЇ ді? 
12С ос "У 
т оба), вай ЕФ 
г дк да 
де Р),  - коефіцієнт | дифузії 
електроліту, а (т) - швидкість появи 
електроліту. 


2.3. Модель біосенсора, яка вико- 
ристовує кінетику Міхаеліса-Ментена 

У роботі |19| розглядається модель 
біосенсора, яка включає в себе три області: 
ферментний шар (ферментна мембрана), де 
відбувається ферментативна реакція; шар 
транспортування  масопередачі, шляхом 
дифузії обох сполук (субстрату 5 1 про- 
дукту Р), через напівпроникну мембрану; 
конвективна область, де концентрація ана- 
літу підтримується постійною. 

Припускаючи симетричну геометрію 
обох мембран та однорідний розподіл 
іммобілізованого ферменту, в фермент- 
ному шарі доцільно використати двоком- 
партментну математичну модель. Дина- 
міку концентрацій субстрату 5, а також 
продукту Р у ферментному шарі можна 
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описати системою рівнянь реакції-дифузії 
(:»0), 


0 
2 

Ор, зар д пе Отак З е , х є (0,1,), 

ОЇ "0 Ки, 


де Х та І - зміщення і час; 5, (3,1) і 


боді і 
дх' Ки 55, 


(9) 


р.(х,І) - концентрації субстрату та 
продукту в ферментному шарі; 4, 1 4р - 
коефіцієнти дифузії; о, - Максимальна 


пах 
ферментативна швидкість (0... З Код» бо - 
концентрація ензиму); К, - константа 
Міхаєліса (К,, -(К | К,)/К, ). 

За межами ферментної мембрани 
відбувається лише транспортування маси 
за рахунок дифузії обох сполук (с»0), 


са, би 

ої" ПОРОДИ 

я й 10) 
коні Зо БА ДЕ М 


де 5, (51), Р (хг) - концентрації 
субстрату та продукту реакції; 4. , 4; - 


коефіцієнти дифузії. 

2.4. Модель біосенсора на основі 
рівнянь реакції-дифузії 

Зв'язування каталітичної реакції 
каталази та пероксидази у ферментатив- 
ному шарі з дифузією в одновимірному 
просторі, описане законом Фіка, приводить 
до наступних рівнянь реакції-дифузії, що 
описують роботу  біосенсора |20| у 
ферментативному шарі (0 «учйзоТа о): 


ОЗ. дз, 
ді гі, вен й 
11 
КаКооваб оба ом 
ПУ ЛИ Кава 
2 
08 ої д'8зг | Коївоавав оба 12 
до од к га а 
Ї Х 2151 Ж Каб е 
др, д'р, ке 
паро -ечочаочіть, 13 
ді зако 2 3) 
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2 
Ор». еВ 0 Ре ц; Коба еба 
ді те 


, 14) 
ок ЛИ ям Коба 


де х та ї означають зміщення і час; 


по); 5».(51) - молярні 
концентрації субстратів 5, 5»; 
Ре (о і ) ; Ре (5 Ї ) - молярні 


концентрації продуктів реакції Р, Р); 


а, - товщина ферментного шару; 


е 


0, З р, і В. 2 РН - дифузні 


коефіцієнти, які вважаються сталими; 
Ку Ко» Ко» - сталі швидкості реакції. 
За | межами ферментного шару 
відбувається тільки масове перенесення 
шляхом дифузії субстратів та продуктів 


реакції. Рівняння керування для 
зовнішнього дифузійного шару 
(120, а «х«а,ча,) представлені 


припущенням, що зовнішний масовий 
транспорт підпорядковується обмеженому 
дифузійному режиму: 


дес д'є, 
вно ря , соб, 5» Р, р. 15) 
де 


за 5» за 5» роб б» | Ра (6 1) 

означають молярні концентрації субстратів 
та продуктів | у з дифузійному шарі 
товщиною 4,; Р, 0, 0,0 - 


а? 7 Ра 
дифузні коефіцієнти. 
Дифузійний шар розглядається як 
дифузійний шар Нернста |21). За підходом 
Нернста шар товщиною а, залишається 


незмінним з часом. Зроблено припущення, 
що розчин знаходиться в русі з однорід- 
ною концентрацією у дифузійному шарі. 

3. Модель імуносенсора у вигляді 
решітчастих деференціальних рівнянь 

У о моделях решітчастих  диферен- 
ціальних рівнянь просторова структура має 
дискретний характер. Решітчаста динаміка 
широко використана в задачах |22, 23, 
оскільки середовища, в яких популяційні 
види живуть, можуть бути дискретними, а 
не неперервними. 
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341. Решітчаста модель взаємодії 
антиген-антитіло для двовимірного ма- 
сиву біопікселів 

Нехай У, (5) - 


антигенів, Б РО; - концентрація антитіл у 


концентрація 


біопікселі (і, )), і, /-1,М. 

Модель грунтується | на | таких 
біологічних припущеннях для довільного 
біопікселя (і, )). 

1. Маємо деяку константу народжуваності, 
В.» 0 для популяції антигенів. 

2. Антигени нейтралізуються антитілами з 
деякою ймовірнісною швидкістю у»б. 

3. Популяція антигенів прагне до деякої 
межі насичення зі швидкістю 0,20. 


4. Ми маємо деяку дифузію антигенів з 
чотирьох сусідніх пікселів (1-1,/), 
(1-1,)), (,)-1), (,)-1) зі швидкістю 


дифузії ДА", де Рр»0 і Л»б0 є 
відстань між пікселями. 

5. Ми маємо деяку сталу смертності 
антитіл ш, 20. 

6. У результаті імунної відповіді ми 
збільшуємо щільність антитіл з 
ймовірнісною швидкістю пу. 

7. Популяція антитіл прагне до деякого 
рівня насичення зі швидкістю 0,20. 


8. Імунна відповідь з'являється з деякою 
сталою затримкою у часі 7»0. 

На основі цього ми розглядаємо дуже 
просту конструкцію моделі  антиген- 
антитіло із запізненням для двовимірного 
масиву біопікселів, яка грунтується на 
добре відомій моделі Марчука |11|, |24, 
25), і використовує просторовий оператор 


дл 


5 , запропонований у (121: 


аУ, (г) 
по Ро 

ай 
-бУ оту, чн 5 У (16) 
4Б, (1) а Еге Ї 
в Лот) 


-бв ФВ, Б, (1) 
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Модель (16) задана початковими 
функціями (17): 
У, (9 а УЮ 20, 


Е (ФВ зБУС20, гє|-г,0), (17) 
У, (0), Е,(0)» 0. 


Побудова моделі грунтувалася на ряді 
біологічних припущень щодо взаємодії 
колоній антигенів та антитіл (26-28). 

Висновки та перспективи подаль- 
ших досліджень 

У роботі розглянуто підходи до роз- 
робки математичних моделей біосенсорів 
та імуносенсорів у вигляді диференціаль- 
них рівнянь у частинних похідних, які б 
враховували біологічні припущення щодо 
основних компонент пристроїв - симет- 
ричну геометрію біосенсорів і однорідний 
розподіл іммобілізованого ферменту та 
реакції зв'язування у ферментному шарі, 
кінетичні властивості бактеріальних фер- 
ментативних реакцій, кінетику Міхаєліса- 
Ментена, колонії антигенів 1 антитіл, 
локалізованих у пікселях, дифузії колонії 
антигенів між пікселями. 

У роботі запропоновано модель 
імуносенсора, яка грунтується на системі 
решітчастих диференціальних рівнянь із 
запізненням. Побудова моделі грунтува- 
лася на ряді біологічних припущень щодо 
взаємодії колоній антигенів та антитіл, а 
також дифузії антигенів. Для опису 
дискретних у просторі колоній, локалі- 
зованих у відповідних пікселях, викорис- 
тано апарат решітчастих диференціальних 
рівнянь. 

У подальших дослідженнях необхід- 
но провести дослідження стійкості моде- 
лей біосенсорів та імуносенсорів на основі 
частинних похідних та решітчастих дифе- 
ренціальних рівнянь із запізненнням. 
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М.У. Когоаїіїі 

АБоші арргоасре5 оп таїпетагбісаї 
5ітшабйоп ої Біо5еп5ог апа іттиповбепяог 
дупатіс 5у5кетя 

Те рарег сопзідег5 Ше арргоасі (0 Ше 
таФетайса! тодеі5 ої біо5еп5ог5 апа 
іпотипобеп5ог5 (паї | (аке | іпіо | ассошпі 
ріоіовіса|  аззигарйоп5 | абоці Ше  таїп 
сотропепі5 ої Ше демісе: Ше 5упітеїгіс 
збеотеїту ої Ше  Біозеп5ог  апд Ше 
ротохепеоця ді5іібийоп ої Фе іптобійлед 
епгуте апа Фе біпаше геасйоп іп Фе епгуте 
Іауег, Ше Кіпейс ргорегііе5 ої Ше Басіепіа! 
епгутайс геасйоп5, Ше Місреїсі8-Мепіеп 
Кіпейся, Ше соіопіе5 апйєєпя апа апібоадїе5 
ІосаПйлед їп ріхеїз, Фійизіоп ої соЇопіе5 ої 
апйяеп5 Бебмееп ріхеїз, апа офегз. 

І Фе могк Ше тоайе! ої іпатипобепзог 
15 ргорозед, уупісі 15 Ба5ед оп Фе зузіет ої 
Іашйсе дїНегеппа! едиайоп5 улії 4сіау. ТРре 
соп5ігисйоп ої Фе тодеі 15 ба5ед оп а питрег 
ої ріоіовіса! аз5итарйопя абоці Фе іпіегасйоп 
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ої соіопіе5 ої апйгеп5 апа апибоаціез, а5 меії! 

а5 Ше ФіНизіоп ої апієеп5. 

Матеїу, Феу аге Ше ГоПом/пе: 

- оуе фБауе зоте соп5апі Бігіргаїс ЙЇог 
апйєеп роршацоп; апйєеп5 аге 
пецітайгед бу апйбодіе5 мулі 50те 
ргобабішу гаге; 

-. апйреп5 аге пецігайгед бу апибоадіе5 мі 
зоте ргобабійсу га/е; 

-  уге рауе 5оте ФфНизіоп ої апібодіе5 гот 
Гоиг пеїєрбогіпє ріхеїз уліб ФіНиз10п гаї; 

- оме Бауе опе дсоп5апі аігібгаїе ої 
апиБоаіе58; 

- 5 а гезції ої шатипе ге5зроп5е ме Баує 
іпсгеазе ої деп5у ої апйбодіе5 мулі 
ргобабішу гак; апабоду рориіайоп (епав 
го 5оте сагтуїпе сарасіїу мій 50те гаї; 

-  шопипе гезроп5е арреаг5 мій 5о0те 
сопуапі те 4еау. 

Тре тобі іпрогіапі Шіпє 15 (0 (аКе іпіо 
ассойпі 5расіаПу із5сгеїе сПагасіег ої Ше 
тоадеі. Оп Фе оФег Папа ме ц5е Фе ргедаюог- 
ргеу тоде! Їог Ше дезсгірйоп ої іттипе 
гезроп5е їп еасп ріхе!. Тре те дсіау ої 
йптипе  гезроп5е 15  езсгібед бу Фе 
аррагация ої деіЇау дійегеппа! едиайопзя. Бог 
Фе Фезсгірйоп ої дївсгеїе зрасез іп Ше 5расе 
ої Фе соіопіе5, ІосаП7ед їп Ше сопезвропдпе 
ріхеї8, Ше аррагаїш5 ої Іакісе діНегепиаї 
едиайоп5 15 | ицзед. ГІакісе  4іНегеппаї 
едиайоп5 агі5е іп таапу аррПеа зибіесівє, зисп 
а5  сретісаї геасйоп, ітаєбе ргосез5іп5, 
тагегіа! 5сіепсе, апа Ббіоїову. п (Пе піоадеізя 
ої Іакісе діНегепца! едиайопя, Ше врайа! 
5ігосгиге Баз а дізсгеїе срагастег, апа Іайісе 
дупатіс5 Паме гесепйу Бееп ехіеп5їуеїу єс 
іо тиоде! Біоіовіса| ргобіетяв зіпсе Ше 
епуїгоптепі їп у/рісп Ше 5ресіе5 роршіайоп 
Пуе5 тау Бе а4ізсгеїе Би пої сопіїпиоц5. 
Усі арргоасі 15 арргоргіаїе ог ошг ргобіет. 
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